A Novel Approach for Evaluating the Influence of Texture Intensities on the First Magnetization Curve and Hysteresis Loss in Fe–Si Alloys

This paper examines the relationship between the magnetization behavior and crystal lattice orientations of Fe–Si alloys intended for magnetic applications. A novel approach is introduced to assess anisotropy of the magnetic losses and first magnetization curves. This method links the magnetocrystalline anisotropy energy of single crystal structures to the textures of polycrystalline materials through a vectorial space description of the crystal unit cell, incorporating vectors for external applied field and saturation magnetization. This study provides a preliminary understanding of how texture influences magnetic loss rates and the first magnetization curves. Experimental results from Electron Back-Scattered Diffraction (EBSD) and Single-Sheet Tests (SSTs), combined with energy considerations and mathematical modeling, reveal the following key findings: (i) a higher density of cubic texture components, whether aligned or rotated relative to the rolling direction, decreases magnetic anisotropy, suggesting that optimizing cubic texture can enhance material performance; (ii) at high magnetic fields, there is no straightforward correlation between energy losses and polarization; and (iii) magnetization rates significantly impact magnetization loss rates, highlighting the importance of considering these rates in optimizing Fe–Si sheet manufacturing processes. These findings offer valuable insights for improving the manufacturing and performance of Fe–Si sheets, emphasizing the need for further exploration of texture effects on magnetic behavior.


Introduction
Magnetic materials are extensively used in commercial applications, including transformers, electrical motors, solenoids, and generators.Among the soft magnetic materials, Non-Oriented Grain (NGO) Fe-Si alloys are the most widely used due to their cost-effectiveness and availability [1][2][3][4][5][6][7].Recently, the demand for higher magnetic performance in automotive electric motors has driven research towards reducing the magnetic energy loss in these steels when subjected to an external magnetic field.
It Is well known that microstructural features such as mean grain size, defects, crystallographic texture, and precipitates influence the magnetic domain movement and, consequently, the magnetic behavior of the material [1][2][3][4][5][6][7].Researchers are thus intensively studying these effects, particularly how to improve the microstructure of Fe-Si alloys during production.
Several studies have highlighted that in Fe-Si steels, crystallographic textures significantly impact the permeability of the first magnetization curve and the area of the hysteresis loop.Specifically, a higher number of crystals with <111> directions aligned with the external applied field results in lower magnetization at a given field intensity and reduced maximum permeability.Conversely, the <100> directions have the opposite effect [1,3,5,7,8].Research on hysteresis loss has considered various aspects of material textures, examining losses at specific values of external fields or magnetic polarization.Key parameters include the following: a generic energy texture parameter, which correlates linearly with loss [9]; average magnetocrystalline energy [10]; texture parameter A θ , dependent on the intensities of relevant magnetic texture components [11]; texture parameter A α , related to the rolling direction and magnetization vector [12]; anisotropy parameter ε [13]; and weighted A θ parameter [14].
This paper introduces a novel method for studying the influence of crystalline anisotropy on the first magnetization curve and total magnetic energy loss.Specifically, the total magnetic energy loss curve is derived from the areas of the hysteresis loops, with the peak polarization values at given external magnetic fields defining the first magnetization curve.It is important to note that magnetic behavior varies between the rolling and transverse directions.
The proposed method includes the following: (i) saturation magnetization; (ii) external applied field vectors; and (iii) energies associated with these vectors within the vectorial space of the crystal unit cell.
This method is implemented through the following steps: (i) Electron Back-Scattered Diffraction (EBSD) analyses and Single-Sheet Tests (SSTs); (ii) EBSD data processing using an appropriate toolbox; (iii) energy considerations related to the magnetization process; and (iv) mathematical modeling of energy loss.
Figure 1 presents a conceptual map of the work.
Research on hysteresis loss has considered various aspects of material textures, examining losses at specific values of external fields or magnetic polarization.Key parameters include the following: a generic energy texture parameter, which correlates linearly with loss [9]; average magnetocrystalline energy [10]; texture parameter  , dependent on the intensities of relevant magnetic texture components [11]; texture parameter  , related to the rolling direction and magnetization vector [12]; anisotropy parameter  [13]; and weighted  parameter [14].
This paper introduces a novel method for studying the influence of crystalline anisotropy on the first magnetization curve and total magnetic energy loss.Specifically, the total magnetic energy loss curve is derived from the areas of the hysteresis loops, with the peak polarization values at given external magnetic fields defining the first magnetization curve.It is important to note that magnetic behavior varies between the rolling and transverse directions.
The proposed method includes the following: (i) saturation magnetization; (ii) external applied field vectors; and (iii) energies associated with these vectors within the vectorial space of the crystal unit cell.
This method is implemented through the following steps: (i) Electron Back-Scattered Diffraction (EBSD) analyses and Single-Sheet Tests (SSTs); (ii) EBSD data processing using an appropriate toolbox; (iii) energy considerations related to the magnetization process; and (iv) mathematical modeling of energy loss.
Figure 1 presents a conceptual map of the work.This approach could enable the development of more accurate models for predicting the magnetization curve and magnetic energy losses compared to existing methods.Furthermore, it will provide preliminary data to help establish guidelines for optimizing rolling processes and heat treatments in industrial plants.This approach could enable the development of more accurate models for predicting the magnetization curve and magnetic energy losses compared to existing methods.Furthermore, it will provide preliminary data to help establish guidelines for optimizing rolling processes and heat treatments in industrial plants.

Theoretical Bases
In the initial part of the first magnetization curve up to the "knee region", as shown in Figure 2, the magnetization process is primarily governed by the motion of domain walls.As the curve progresses through the "knee region", the saturation magnetization vectors of the newly formed domains begin to rotate, aligning with the direction of the externally applied magnetic field.Eventually, at high external magnetic field values, the domain configuration reaches a state where all saturation magnetization vectors are parallel to the applied field direction.

Theoretical Bases
In the initial part of the first magnetization curve up to the "knee region", as shown in Figure 2, the magnetization process is primarily governed by the motion of domain walls.As the curve progresses through the "knee region," the saturation magnetization vectors of the newly formed domains begin to rotate, aligning with the direction of the externally applied magnetic field.Eventually, at high external magnetic field values, the domain configuration reaches a state where all saturation magnetization vectors are parallel to the applied field direction.For a single bcc crystal structure material subjected to an external magnetic field applied along one of the crystal directions <100>, <110>, and <111>, the magnetization value at the initial part of the "knee region" varies as follows: it is highest for <100>, lowest for <111>, and intermediate for <110>.However, the external magnetic field required to align all saturation magnetization vectors parallel to the applied field is lowest for <100>, highest for <110>, and intermediate for <111>.
The mean magnetization of a material is related to the variation in the Landau free energy.Specifically, under constant entropy and temperature conditions, it is given by the Helmholtz free energy, which reflects changes in the system's internal energy.
In an ideal bcc crystal, the potential energy of the saturation magnetization vector  , which drives the magnetization process up to the "knee region," is given by: where   is the saturation magnetization vector;  is the external applied magnetic field vector; and  is the vacuum permeability; the symbol denotes the scalar product between vectors.However, the magnetization process up to saturation is also influenced by the magnetocrystalline anisotropy energy  .For a single bcc crystal structure material subjected to an external magnetic field applied along one of the crystal directions <100>, <110>, and <111>, the magnetization value at the initial part of the "knee region" varies as follows: it is highest for <100>, lowest for <111>, and intermediate for <110>.However, the external magnetic field required to align all saturation magnetization vectors parallel to the applied field is lowest for <100>, highest for <110>, and intermediate for <111>.
The mean magnetization of a material is related to the variation in the Landau free energy.Specifically, under constant entropy and temperature conditions, it is given by the Helmholtz free energy, which reflects changes in the system's internal energy.
In an ideal bcc crystal, the potential energy of the saturation magnetization vector E p , which drives the magnetization process up to the "knee region," is given by: where M s is the saturation magnetization vector; H is the external applied magnetic field vector; and µ 0 is the vacuum permeability; the symbol denotes the scalar product between vectors.However, the magnetization process up to saturation is also influenced by the magnetocrystalline anisotropy energy E a .
Here, α i represents the direction cosine, which is the angle between the saturation magnetization vector and a specific crystal direction during its rotation toward the direction of the external applied field, where it reaches the minimum value of E a .Specifically, the magnetocrystalline anisotropy energy E a follows the order E a,<100> < E a,<110> < E a,<111> when the saturation magnetization vectors are fully aligned with the external applied field [1,2,5,7,8].
Moreover, in materials like iron or Fe-Si alloys, the term K 2 could be neglected since [1,8,15].
The total energy of the ideal material during the magnetization process up to the "knee region" (E t,1 ), neglecting other energy contributions like the shape one, is as follows: While over the "knee region", considering also the anisotropy magnetocrystalline energy, the total energy (E t,2 ) is calculated as follows: In real polycrystalline materials, the magnetization process is more complex due to factors such as the frequency of the external magnetic field and the specific crystal orientations of each grain.Additionally, materials often contain various defects-such as non-metallic inclusions, voids, and second phases-that act as pinning sites for domain wall motion.Despite these complexities, this study focuses on crystal orientations.Consequently, factors such as frequency, shape anisotropy, and inclusions are considered constant, as they are consistent across the specimens studied.
In polycrystals, lattice orientations can be described using Euler angles, which represent the rotation of the reference frame fixed to the unit crystal cell relative to the reference frame fixed to the specimen.The most comprehensive information about the distribution of crystal orientations in the material is provided by orientation distribution function (ODF) maps [15][16][17][18][19].
The ODFs for each orientation g are given by the orientation function distribution density f (g).
where dV is the volume of the crystal with orientation g; V is the total volume of the crystal; and dg is the orientation variation from g [15].It is also possible to plot the ODF along a particular set of directions, obtaining a so-called fiber plot.
To assess the relationships between energy contributions and texture analyses, two reference frames are considered: first, the sheet reference frame (RD, TD, and ND), as shown in Figure 7 (d2); second, the specimen reference frame (X s , Y s , Z s ), as shown in Figure 7 (d3).The reference frame axes are aligned such that RD ∥ X s , ND ∥ Y s , and TD ∥ Z s .
The most important rolled fibers for bcc crystal materials are α, with < 110 >∥ X s ∥ RD, and γ, with < 111 >∥ Z s ∥ TD.Specifically, the orientations are defined by a set of Euler angles that describe the rotation from the initial orientation g i to the final orientation g f .For α and γ fibers the set of orientations are (Table 1): In recrystallized Fe-Si alloys, the commonly considered texture components include the Goss component {110}<001> and the cube component {100}<001> [16,17,19,20].How-ever, this study will focus on analyzing the set of cubic orientations, also known as the "fiber Cube" (Table 2).
Cube : [001] ∥ Z s ∥ TD, considering all the rotated cubes with respect to the X s ∥ RD and considering the magnetic behavior along the [100] and [010] axes.
This assumption enables the calculation of E t,2 along each principal crystal axis while varying only one Euler angle.Table 3 provides the anisotropy energy E a and potential energy E p along the rolling direction (RD): Table 3. Anisotropy (E a ) and potential energies (E p ) along RD.
In this context, β 1 is the angle between the external applied field H and the saturation magnetization M s along one of the easy axes (<100>), with 0 ≤ β 1 ≤ 45 • , accounting for the symmetry of M s along the crystal axes; θ 1 is the angle between M s and the <100> directions, with 0 ≤ θ 1 ≤ 45 • ; θ 2 is the angle between M s and the <100> directions, with 0 ≤ θ 2 ≤ 45 • ; and δ 1 is the misorientation angle between M s and H, given that H belongs to the {111} planes.
The definitions of the angles θ 1 and θ 2 imply the following relationships: Thus, for the anisotropy energy: E a,γ < E a,Cube < E a,α .Therefore, the α fiber is the one that generates the highest value of magnetocrystalline anisotropy energy opposing the alignment with the external magnetic field.
The potential energy is influenced by the Thus, it follows that: (i) the fewer the <100> directions that are misoriented with respect to the rolling direction (RD), the higher the potential energy of the magnetic moments and the lower the anisotropy field.As a result, the magnetization value at the beginning of the "knee region" is higher, and saturation magnetization can be achieved with a lower external magnetic field.(ii) The E a contribution of γ fiber can be neglected because K 0 ∼ = 0 and only E p ̸ = 0. (iii) E a, α provides the highest energetic contribution.The higher its value, the lower the magnetization value at high external magnetic field intensities H.
Table 4 reports anisotropy E a and potential E p along TD.
Table 4. Anisotropy (E a ) and potential (E p ) energies along TD.
Where β 2 is the angle between H and M s along one of the easy axes, with 0 ≤ β 2 ≤ 45 • due to the symmetry of M s along the crystal axes; θ 3 is the angle between M s and the <100> directions, with 0 ≤ θ 3 ≤ β 2 ; and δ 2 is the misorientation angle between M s and H, where H is parallel to <111>.
The definition of the angles θ 3 implies that E a,Cube < E a,α < E a,γ .Thus, the γ fiber generates the highest magnetocrystalline anisotropy energy, opposing the alignment with the external magnetic field.
Regarding the potential energy: if 0 ≤ β 2 < 45 • , then E p,Cube ≥ E p,α .Thus, the γ produces the highest magnetocrystalline anisotropy energy opposing the external applied field.This results in a reduced magnetization value at the beginning of the "knee region", and the saturation magnetization can only be achieved at higher external magnetic field values compared to those along the rolling direction (RD).Additionally, the E a contribution of the "Cube" fiber can be neglected because K 0 ∼ = 0 and E p, Cube = max E p .
Ultimately, because E a, γ provides the highest anisotropy contribution among the fibers, a higher value of E a, γ results in a lower magnetization value at high external magnetic field intensities H.
From these observations, anisotropy energy contributions can be assessed for an external magnetic field applied along specific crystal directions or within particular crystal planes.This leads to a simplification in studying the effects of anisotropy energies.
To qualitatively assess the influence of textures on energy contributions, we consider a simplified analysis of the function f (g) along a fiber.
Given that g ∈ R 3 , f : R 3 → R , and since two angles remain constant while one angle varies along the fiber, we can define a function h(θ) where θ is the varying angle along the fiber, with θ ∈ R and h : R → R .
The function h(θ) is not defined but is derived from f (g) along a fiber.
To account for all fiber orientation contributions with one specific direction parallel to the specimen-fixed reference frame, a numerical integration of the function h(θ) − I is performed from the initial varying angle θ i to the final angle θ f : The study of the influence of crystallographic anisotropy on total energy loss can be conducted by describing the phenomenon using a mathematical function, following an identification process [21].This function does not directly represent the physical processes of magnetism but provides a mathematical model useful for a preliminary analysis of how various factors affect the material's behavior.
Given the shape of the total energy loss, a double exponential function is selected as the fitting model: In this study, the rates of variation of total energy loss up to the "knee region" and in the rotating regions were chosen since the function could be linearized in two main regions.
The first region is around H = 0 A m , with a McLaurin series until the first degree: The coefficient of the linearized function, which represents the rate of variation in the first region, is: The second region corresponds to the high value of H where the function can be approximated with a line with starting point P 1 = H ≡ H re f ; P s, f it H re f and ending at point P 2 = H ≡ H end = 10000 A m ; P s, f it (H end ) : The linear approximating function is as follows: The coefficient of the linearized function, which represents the rate of variation in the second region, is as follows: m lin,2 = p 1 (11)

Materials and Methods
The materials studied were three rolled Fe-Si alloys of grade M350-50A, in accordance with EN 10106.Their chemical compositions were analyzed using an Optical Emission Spectrometer (OES ARL 346 (Thermo Fisher Scientific, 168 Third Avenue, Waltham, MA USA 02451)).The acronyms and compositions of the analyzed sheets are detailed in Table 5.The key alloying elements are C, Si, Mn, and Al, as they significantly impact the magnetization process, saturation polarization, energy loss, and magnetocrystalline anisotropy.Additionally, S and N influence the total energy loss [1,2,6,8,10].Magnetic characterization was performed using a Brockhaus Messtechnik Single-Sheet Tester (SST) (BROCKHAUS MEASUREMENTS, Dr. Brockhaus Messtechnik GmbH & Co. KG, Lüdenscheid, Germany) model MPG100 D DC/AC.The device features an external magnetic field frequency range from 3 Hz to 10 kHz, with a maximum polarization of 2 T and measurement repeatability within ±2%.
The external magnetic field was applied along the rolling direction (RD) and transverse direction (TD) of the sheets at a frequency of 50 Hz.
Before investigation, samples were embedded in conductive resin and polished according to standard metallographic procedures, achieving a finish with 50 nm colloidal silica.Table 6 summarizes the differences in magnetic energy loss along RD and TD at  .L1 exhibits a negative difference, indicating higher energy loss along TD compared to RD, while L3 has the highest difference between RD and TD.The energy loss curves are derived from the value areas of hysteresis loops, whose magnetic polarization tips and corresponding external magnetic field intensities generate the first magnetization curve (Figure 5).

Magnetic Characterization
Figure 4 presents the experimental energy losses for each material along the respective RD and TD.The curves show that P s is higher at H end only for L1 along TD than RD.Table 6 summarizes the differences in magnetic energy loss along RD and TD at  .L1 exhibits a negative difference, indicating higher energy loss along TD compared to RD, while L3 has the highest difference between RD and TD.Table 6 summarizes the differences in magnetic energy loss along RD and TD at H end .L1 exhibits a negative difference, indicating higher energy loss along TD compared to RD, while L3 has the highest difference between RD and TD.Table 7 provides the polarization J[T] at  .The data show that ( ) > ( ) > ( ) along RD while ( ) > ( ) > ( ) along TD.
Table 7. Polarization values at the maximum value of the externally applied field along RD and TD.Table 8 reports the differences in polarization Δ along RD and TD at  =  = 700

𝐽[𝑇]
, considered the initial point of the "knee region."It shows that polarization along RD is consistently higher than along TD at  .Table 7 provides the polarization J[T] at H end .The data show that J(L 3 ) > J(L 2 ) > J(L 1 ) along RD while J(L 1 ) > J(L 2 ) > J(L 3 ) along TD.
Table 7. Polarization values at the maximum value of the externally applied field along RD and TD.

J[T]
RD TD L1 Figure 6 presents the experimental first magnetization curves for each material along RD and TD.The graphs show that for all the sheets J[T] at H end is higher along RD compared to TD.Table 7 provides the polarization J[T] at  .The data show that ( ) > ( ) > ( ) along RD while ( ) > ( ) > ( ) along TD.
Table 7. Polarization values at the maximum value of the externally applied field along RD and TD.Table 8 reports the differences in polarization Δ along RD and TD at  =  = 700

𝐽[𝑇]
, considered the initial point of the "knee region."It shows that polarization along RD is consistently higher than along TD at  .Table 8 reports the differences in polarization ∆J along RD and TD at H = H k = 700 A m , considered the initial point of the "knee region".It shows that polarization along RD is consistently higher than along TD at H k .Table 9 presents the differences in polarization between all sheets along RD and TD at H k .The results indicate that L2 shows the greatest differences in magnetization compared to the other sheets.The fitting parameters for Function (7) were obtained (Table 11) by solving a nonlinear least-squares problem [21,24] using various types of functions and employing a Trust Region Algorithm [24] via the MATLAB R2024a Curve Fitting Toolbox.The goodness of fit for the experimental data was evaluated using the sum of squared errors (SSE) [24] and the root mean squared errors (RMSE) [25] to determine if the Trust Region Algorithm successfully converged to the optimal minimum.Additionally, the coefficient of determination R 2 was calculated to assess the correlation between the experimental data and the fitting curve [26], with particular emphasis on the adjusted R 2 .

Data Analyses
This adjusted R 2 accounts for the number of fitting parameters and helps identify potential overfitting [27].
The low values of SSE and RMSE, along with very high values of R 2 (Table 12), confirm a good fit between the experimental and calculated data.Thus, Function (7) accurately describes the magnetic energy loss behavior of the material along the rolling direction (RD) and transverse direction (TD) and is suitable for studying the material's magnetic anisotropy as a function of the intensity of the applied external field.Based on the above, Function (7) can be reformulated as follows: The linear coefficient values of Function (8) are reported in Table 13 for all the energy loss fitted curves, indicating that  >  for all the sheets.(10) accurately fits the data provided by Function (7) for  ≥  , while the  values confirm a strong correlation between the experimental and fitted data.Based on the above, Function (7) can be reformulated as follows: The linear coefficient values of Function (8) are reported in Table 13 for all the energy loss fitted curves, indicating that  >  for all the sheets.(10) accurately fits the data provided by Function (7) for  ≥  , while the  values confirm a strong correlation between the experimental and fitted data.Based on the above, Function (7) can be reformulated as follows: The linear coefficient values of Function ( 8) are reported in Table 13 for all the energy loss fitted curves, indicating that m RD > m TD for all the sheets.The procedure for finding the best line approximation is the same as applied to Function (7).The values of linear coefficients for the fitting linear Function (10) are reported in Table 15, showing that m RD < m TD for all the sheets.

Discussion
Based on the results of the first magnetization curves along the rolling direction (RD) and transverse direction (TD) (Figures 5 and 6), and considering Equations ( 1) and (3) along with Tables 3, 4 and 8-10, the following conclusions can be drawn for H ≤ H k where domain wall motion is the predominant mechanism: where E p,Cube = −M s H; E p,α = −M s Hcos(45 • ); E p,γ = −M s Hcos(54.7 • ). Thus, At H k , considering the integral of textures along RD and TD: Hence, However, ), the significant difference I Cube (L 3 ) ≫ I Cube (L 1 ) explains the results.These observations are supported by polarization differences along RD and TD for H = H k .
Figure 5 points out that the first magnetization curve of L2 consistently shows higher polarization values than L1 and L3, with L3 being higher than L1, and L1 having the lowest polarization values among the curves.Table 9 shows these differences for H = H k , and the texture integral values from Table 10 confirm the results from the previous equations.
Figure 6 highlights that for H = H k the polarization values are consistently lower along TD than RD.Table 8 illustrates these differences for H = H k , and the texture integral values from Table 10 further confirm the results from the previous equations.
For intensities of external applied field lower than H k for each sheet, denoted as H low , the following expression is used: These observations are supported by polarization differences along RD and TD for H < H k .
From Figure 6, it is evident that from very low values of external applied field intensity, the polarization values are consistently lower along TD than RD.The texture integral values from Table 10 verify the results from the previous equations.
These observations demonstrate that higher Cube texture integral values correspond to higher magnetization at the onset of the "knee region", indicating that less energy is required for magnetization with this texture and that magnetocrystalline anisotropy also plays a role at relatively low magnetic fields.
For H → H end , considering Equation (4), Tables 3, 4, 7 and 10, along with Figures 5 and 6, where E p is approximately constant for all sheets, the following expression is used: Along RD: The following is given: Thus, These results are confirmed by polarization values at maximum external fields along RD. Figure 5a highlights that for H → H end , the polarization values are the highest for L3, lower for L2, and the lowest for L1.Table 7 shows these differences for H → H end , and the texture integral values from Table 10 prove the results from the previous equations, confirming that lower Cube and α texture integral values correspond to higher magnetization at maximum magnetic field intensity.
Along TD: Given that E a,Cube ∼ = 0, E a,TD = E a,α + E a,γ . Since Thus, E a,TD (L 3 ) > E a,TD (L 2 ) > E a,TD (L 1 ). Therefore, These results are confirmed by polarization values at maximum external fields along TD. Figure 5b points out that for H → H end along TD, the polarization values are the highest for L1, lower for L2, and the lowest for L3.Table 7 shows these differences for H → H end , and the texture integral values from Table 10 support the results from the previous equations, confirming that lower Cube and γ texture integral values correspond to higher magnetization at the maximum field intensity.
Hence, anisotropy energy is lower at high magnetic field intensities, and the material exhibits less isotropic behavior.
According to [1,5,7], increasing silicon content in Fe-Si alloys generally increases electrical resistivity and should decrease energy loss.However, experimental data show the highest energy loss for L1 (1.879% Si), the lowest for L2 (1.718% Si), and intermediate for L3 (1.971% Si) at H end along RD and TD (Figure 1).This discrepancy suggests that the magnetic behavior of Fe-Si alloys should also consider the effects of N, S, Al, and Mn (Table 5).Specifically, elements like AlN can increase magnetic energy loss.Given the low and nearly constant S content, the presence of AlN is significant: L1, with the highest Al and N content (0.403% Al, 0.0086% N), exhibits the worst magnetic behavior; L3, with low Al and intermediate N (0.353% Al, 0.0042% N), shows intermediate behavior; and L2, with intermediate Al and the lowest N (0.397% Al, 0.003% N), demonstrates the best magnetic behavior.
Examining the initial part of the loss curve (Figures 8 and 9) and using Equation (1), Equation (3), and Equation (8), the following can be derived: P s, f it ∝ J(H).Therefore, ∆P s, f it ∝ ∆J(H) and lim The following is given: Thus, µ ′ i ∝ m lin,1 , with µ ′ i being the initial differential permeability.Assuming homogeneous defect distribution and nearly constant electrical resistivity, at H → 0 , the following can be derived: Thus, E a ̸ = 0.
From Figure 8, it is evident that along RD and TD for H ≤ H k the magnetic loss curve of L2 consistently shows higher magnetic loss values than L1 and L3, with L3 being higher than L1, and L1 having the lowest magnetic loss values among the curves.By linearizing the first magnetization curves (Figures 4 and 5) with an equation similar to Equation ( 8) from H → 0 up to H k , to avoid instrumental acquisition errors at very low fields, it is observed that higher J(H k ) correlates with a higher loss rate, as indicated by Equation ( 9) and supported by Table 13.Based on linear coefficients (Table 15) and the previous observations of the first magnetization curve and loss rate for H < H k , anisotropy energy significantly influences magnetic energy loss behavior.
From Figure 9, it is evident that for H ≥ H re f the magnetic loss curve of L1 consistently shows higher magnetic loss values along TD than RD compared to other sheets, while the others always have lower values along TD than RD.Table 15 shows these differences for rates of loss for H ≥ H re f , and together with Table 6 and the texture integral values from Table 10, it is confirmed that texture energies impact L1 more than other sheets for H ≥ H re f along RD and TD.
The results presented represent a novel method for describing magnetic behavior in the vectorial space of the crystal unit cell, involving saturation magnetization, external applied field vectors, and associated energies for typical fiber texture orientations in rolled materials.This approach offers a new perspective compared to semi-empirical methods.Further research is needed to explore the contributions of different orientations to the magnetization and energy loss of materials.

Conclusions
The aim of this paper is to evaluate the effect of material texture on the magnetization and magnetic energy loss of Fe-Si alloys by considering the definitions of potential and magnetocrystalline anisotropy magnetization energy in magnetic materials.This study integrates the orientation distribution function along each fiber and models energy loss as a function of the externally applied field.
The results highlight the following key points: 1. Impact of texture on energy loss: the texture of the material significantly influences the first magnetization curve and the initial energy loss.According to Table 13, the rate of variation in energy loss between the rolling direction (RD) and the transverse direction (TD) is 2.04% for each specimen.Specifically, along the RD, the energy loss rates are as follows: 46.93% between L1 and L2; 15.51% between L1 and L3; and 19.44% between L2 and L3.Along the TD, the rates are: 32.4% between L1 and L2; 15.79% between L1 and L3; and 19.72% between L2 and L3. 2.
Variation in energy loss with polarization: although the maximum polarizations achieved by L3 and L2 are similar along RD and TD, their energy loss behavior differs significantly with the applied external field.Table 7 shows that the polarization difference between L2 and L3 is 0.17% along RD and 0.67% along TD.The variation in energy loss rates between L2 and L3 is 9.09% along RD and 12.64% along TD.Additionally, the polarization difference between RD and TD for the same specimen is 0.0648 T for L2 and 0.0795 T for L3 (Table 7).The corresponding variation in energy losses between RD and TD is 0.0639 W/kg for L2 and 0.2057 W/kg for L3 (Table 6).This results in a 1.38% variation in energy loss relative to the maximum polarization for L2 and a 61.35% variation for L3, indicating that the first magnetization curve does not directly influence the energy loss behavior.

3.
Modeling energy loss with double exponential curve: The influence of texture on energy loss can be effectively studied using a double exponential curve, which provides insight into the contributions of different textures.

4.
Vectorial space analysis: describing saturation magnetization, external applied field vectors, and the related energies in the vectorial space of the crystal unit cell-based on a typical fiber texture set of orientations for the rolled material-proves useful for analyzing both magnetization and magnetic energy losses.5.
These findings underscore the complex interplay between material texture and magnetic behavior, offering a detailed understanding of how different textures affect magnetization and energy loss in Fe-Si alloys.

Figure 1 .
Figure 1.Conceptual map of the paper's work.

Figure 1 .
Figure 1.Conceptual map of the paper's work.

Figure 2 .
Figure 2. Magnetic polarization J vs. external applied magnetic field H and "knee region" highlighted.

Figure 2 .
Figure 2. Magnetic polarization J vs. external applied magnetic field H and "knee region" highlighted.

Figure 3
Figure3illustrates the energy losses P s W kg of all the materials in function of the intensity of the externally applied H A m field along RD and TD, as evaluated by the SST method.It highlights that at H end : P s (L 1 ) > P s (L 3 ) > P s (L 2 ) in RD and TD.

Figure 3 .
Figure 3. Magnetic characterization experimental energy loss  vs. external applied magnetic field H: energy losses along RD (a) and TD (b).The colors are related to the behavior of: L1-blue; L2orange; and L3-yellow.

Figure 4 Figure 4 .
Figure 4 presents the experimental energy losses for each material along the respective RD and TD.The curves show that  is higher at  only for L1 along TD than RD.

Figure 5
Figure 5 displays the experimental first magnetization curves along RD and TD up to  .

Figure 3 .
Figure 3. Magnetic characterization experimental energy loss P s vs. external applied magnetic field H: energy losses along RD (a) and TD (b).The colors are related to the behavior of: L1-blue; L2-orange; and L3-yellow.

Figure 3 .
Figure 3. Magnetic characterization experimental energy loss  vs. external applied magnetic field H: energy losses along RD (a) and TD (b).The colors are related to the behavior of: L1-blue; L2orange; and L3-yellow.

Figure 4 Figure 4 .
Figure 4 presents the experimental energy losses for each material along the respective RD and TD.The curves show that  is higher at  only for L1 along TD than RD.

Figure 5
Figure 5 displays the experimental first magnetization curves along RD and TD up to  .

Figure 4 .
Figure 4. Magnetic characterization experimental energy loss P s vs. external applied magnetic field H: energy losses along RD (blue) and TD (orange).The figures are related to: (a) L1; (b) L2; and (c) L3.

Figure 5
Figure 5 displays the experimental first magnetization curves along RD and TD up to H end .

Figure 5 .
Figure 5. Magnetic characterization polarization J vs. external applied magnetic field H: first magnetization curves along RD (a) and TD (b).The colors are related to the behavior of: L1-blue; L2orange; and L3-yellow.

Figure 6 Figure 6 .
Figure 6 presents the experimental first magnetization curves for each material along RD and TD.The graphs show that for all the sheets [] at  is higher along RD compared to TD.

Figure 5 .
Figure 5. Magnetic characterization polarization J vs. external applied magnetic field H: first magnetization curves along RD (a) and TD (b).The colors are related to the behavior of: L1-blue; L2-orange; and L3-yellow.

Figure 5 .
Figure 5. Magnetic characterization polarization J vs. external applied magnetic field H: first magnetization curves along RD (a) and TD (b).The colors are related to the behavior of: L1-blue; L2orange; and L3-yellow.

Figure 6 Figure 6 .
Figure 6 presents the experimental first magnetization curves for each material along RD and TD.The graphs show that for all the sheets [] at  is higher along RD compared to TD.

Figure 7 Figure 7 .
Figure 7 summarizes the results of EBSD analyses, processed with the MATLAB Toolbox MTEX.It shows that the [111] direction (blue) aligned with the normal direction of the transverse section is predominant in all sheets.Additionally, some grains with [001] (red) and [011] (green) directions are present in L2 and L3. Materials 2024, 17, x FOR PEER REVIEW 11 of 19

Figure 8 Figure 8 .Figure 9 .
Figure 8 reports the fitted curves P s, f it of Figure 3, while Figure 9 reports the fitted curves of Figure 4.

Figure 8 .Figure 8 .Figure 9 .
Figure 8. Fitted curves of experimental data of magnetic energy loss P s, f it vs. external applied magnetic field H along RD (a) and TD (b).The colors are related to the behavior of: L1-blue; L2-orange; and L3-yellow.

Figure 9 .
Figure 9. Fitted curves of experimental data of magnetic energy loss P s, f it vs. external applied magnetic field H along RD (blue) and TD (orange).The figures are related to: (a) L1; (b) L2; and (c) L3.

Table 1 .
Set of Euler angles of fibers α and γ.

Table 2 .
Set of Euler angles of fiber Cube.

Table 5 .
Chemical composition (weight %) of the studied Fe-Si alloys.

Table 6 .
Total energy loss differences between RD and TD for all the sheets.

Table 6 .
Total energy loss differences between RD and TD for all the sheets.

Table 6 .
Total energy loss differences between RD and TD for all the sheets.

Table 8 .
Differences of polarizations between RD and TD for each sheet.

Table 9 .
Differences of polarization along RD and TD between all the sheets.

Table 10
lists the integral values  of each fiber for each sheet.The data show that L2 has the highest  , L3 the highest  , and L1 has the highest  .

Table 10
lists the integral values I of each fiber for each sheet.The data show that L2 has the highest I Cube , L3 the highest I γ , and L1 has the highest I α .

Table 10 .
Values of fiber integrals.

Table 14
provides the SSE, RMSE, and  values for  ≥  , where  is the field strength at which the experimental curves become linear.The SSE and RMSE values indicate that the linear function

Table 14
provides the SSE, RMSE, and  values for  ≥  , where  is the field strength at which the experimental curves become linear.The SSE and RMSE values indicate that the linear function

Table 14
(7))ides the SSE, RMSE, and R 2 values for H ≥ H re f , where H re f is the field strength at which the experimental curves become linear.The SSE and RMSE values indicate that the linear function(10)accurately fits the data provided by Function(7)for H ≥ H re f , while the R 2 values confirm a strong correlation between the experimental and fitted data.

Table 14 .
Values of SSE, RMSE, and R 2 for each sheet for fitted magnetic energy loss with Function (5).
,RD (H low ) = E p,TD (H low )and ∆J > 0 Then, E a,RD (H low ) < E a,TD (H low ) Thus, J RD (H low ) > J TD (H low ).
curation.A.F.: resources and data curation.All authors have read and agreed to the published version of the manuscript.The research was financed by: The European Union-NextGenerationEU (National Sustainable Mobility Center CN00000023, Italian Ministry of University and Research Decree No. 1033-17/06/2022, Spoke 11-Innovative Materials and Lightweighting); Regione Emilia-Romagna-POR FESR 2021-2027 ACTION 1.1.2Call For Strategic Industrial Research Projects Aimed At The Priority Areas Of The Smart Specialisation Strategy DGR 2097/2022"-ID 37993 PG/2023/303085 CUP E27G22000330003 (project MAGMOVE-Sustainable magnetic materials in high-efficiency electric motors for the Green transition).The opinions expressed are those of the authors only and should not be considered as representative of the European Union or the European Commission's official position.Neither the European Union nor the European Commission can be held responsible for them.
Funding:Informed Consent Statement: Not applicable.